Violation of the area law and long range correlations in infinite matrix product states 
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We propose to construct a family of states of one-dimensional spin chains by replacing the finite 
dimensional matrices in matrix product states by infinite dimensional operators constructed from 
bosonic annihilation and creation operators. The resulting states are demonstrated to violate the 
area law and to exhibit long range correlations. In addition, we propose an efflcient way to prepare 
the states experimentally, in which the spins interact sequentially with an ancilla system. 
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I. INTRODUCTION 

Even if the underlying physical model is well-known, it 
is, in general, difficult to simulate quantum many body 
systems on a classical computer. The difficulty stems 
from the fact that the dimension of the involved Hilbert 
space grows exponentially with the size of the system, 
which often renders a treatment in the full Hilbert space 
intractable. This motivates the search for families of 
states that provide a good approximate description of the 
solution to a specific set of problems but are specified in 
terms of a much smaller number of parameters. 

Considering a one dimensional chain of spins, one such 
family of states is the matrix product states (MPS). An 
MPS (see [J for a review) with open boundary conditions 
is a state of the form 



where A^s} are D x D matrices, (u| is a 1 x D vector, \v) 
is a -D X 1 vector, and \si) is the state of the zth spin. 
The number of parameters required to specify the state 
scales linearly in the number of spins N and quadratically 
in the virtual dimension D. Matrix product states are 
particularly suitable to approximate the ground state at 
zero temperature of Hamiltonians with finite range spin- 
spin interactions and an energy gap between the ground 
state and the lowest excited state. This appears because 
these states fuffil the area law 0, 0], which states that 
the von Neumann entropy of the reduced density matrix 
of a subchain is proportional to the number of spins at 
the boundary, i.e., two in one dimension. For an MPS 
the entropy of a subchain is limited by the logarithm of 
D [l| , and it has been shown that all states fulfilling the 
area law can be approximated by MPS using a number 
of parameters that scales polynomially in the number of 
spins 

There are, however, a number of interesting cases, 
where MPSs either fail to approximate the solution or 
provide less accurate results. A particularly interesting 
example of the latter is critical systems, where the en- 
tanglement entrow of a subchain consisting of L spins 
scales as ln(L) WJ^- It has been shown that such systems 



can be approximated by MPSs if the virtual dimension 
grows polynomially in the size of the system fl\ , and the 
multiscale entanglement renormalization ansatz has also 
been used to investigate critical systems @, Q . To ob- 
tain states with a higher entropy than MPSs, it has been 
proposed in [lo| to replace the finite dimensional matri- 
ces A^^^ by infinite dimensional chiral vertex operators 
constructed from a conformal bosonic field [ll|. The re- 
sulting states, which are called infinite matrix product 
states (iMPS) because D — > oo, are able to describe crit- 
ical and noncritical systems on an equal footing. 

The chiral vertex operators contain a product of in- 
finitely many normal ordered exponentials of bosonic an- 
nihilation and creation operators, where the operators 
in different exponentials commute because they corre- 
spond to different modes. Since each of the infinitely 
many bosonic modes constitutes an infinite dimensional 
Hilbert space on its own, it is, in fact, sufficient to in- 
clude just one of them to obtain an iMPS, and the aim 
of the present paper is to investigate the properties of 
the family of states that emerge, when we keep only M 
of the bosonic modes in the chiral vertex operators. We 
find that these states contain even more entropy than the 
critical iMPS in [l^l and exhibit long range correlations. 

A further motivation for investigating this family of 
states is that the particular construction of the states 
naturally suggests a way to prepare the states experi- 
mentally in which each of the spins in the chain interacts 
sequentially with M or A/ + 1 ancilla systems, whereafter 
a conditional measurement is applied to each of the an- 
cillas. The preparation scheme is efficient (for moderate 
M) because the required number of operations per an- 
cilla scales linearly in N and we find that the probability 
of a successful outcome is approximately of order 1/iV for 
each of the conditional measurements. We note that the 
idea of sequential generation is in line with [l^fisj . where 
it is shown that all MPSs can be generated determinis- 
tically by letting a chain of spins initially in a product 
state interact sequentially with an ancilla system with D 
levels, and that all states prepared in this way are in- 
stances of MPSs with virtual dimension D. In [1^, it is 
shown that also continuous matrix product states can be 
generated deterministically by sequential interaction of a 
continuous quantum system with an ancilla. 
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The structure of the article is as follows. In Sec. |TT1 
we provide explicit expressions for the states obtained 
by keeping only M modes in the vertex operators, and 
we compute entropies and correlation functions of these 
states to characterize their properties. In Sec. IIIIl we 
propose and discuss a way to generate the states exper- 
imentally. The proposed experimental implementation 
suggests some natural modifications of the scheme, and 
the influence of these modifications on the properties of 
the produced states is investigated in Sec.|TVl Section IVl 
concludes the paper, and the appendices provide details 
on the derivations of the analytical results stated in the 
main part of the article. 



II. IMPS CONSTRUCTED FROM ONE OR A 
FEW BOSONIC MODES 

A. Wave function 

In general, one can write the state of a one-dimensional 
chain of N spins as 



S2 



sn 



,Sn)\Si, 



(2) 



where Si is summed over all possible states of the ith spin 
and tp{si, . . . ,sn) are complex coefficients. We shall here 
assume Si = ±1. For the family of states investigated in 
the present work, ip{si, . . . , sn) is chosen as 

ip{si, ...,sn) = {Vs,{zi)Vs:,{z2) . ..Vs,^{zn)), (3) 

where 

Vs„ (zn) = exp iSs,i^/a(t)o + s„ E ~/="in^n I 

\ rn—1 ^ / 

X exp (5s„ V^TTo ln(z„) - s„Va ^ ^=a„z„ " (4) 

\ m— 1 ^ / 

and we shall always take N to be even. Here, 5 is 
either or 1, a is a real and positive parameter, z„ 
is a complex parameter, ""Oj and am are operators 
for which the nonvanishing commutation relations are 



hiT^o] = i and [a„,aT 



the expectation value in 



([3]) is evaluated with respect to the state |0) defined by 
ttqIO) = flmlO) = 0, and we note that the resulting state 
is not necessarily normalized. Graphical illustrations of 
the MPS and IMPS constructions are given in Fig. [TJ 

In the limit M — co and (5 = 1, Q is a chiral vertex 
operator of a conformal field theory with central charge 
c = 1 11 1 as investigated in (T^J. In the following, we 
shall refer to the mode on which and ttq act as the 
zero mode and use the term mth mode for the mode on 
which the bosonic annihilation operator am acts. Note 
that 5 = 1 {5 = Q) corresponds to keeping (omitting) 
the zero mode. We consider both cases because the zero 
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FIG. 1. (a) Graphical representation of an MPS as a ten- 
sor network. The tensors Ai] are drawn as rectangles with 
two horizontal legs corresponding to the two virtual indices 
(running from 1 to D) and one vertical leg representing the 
physical index Si. The input and output states are drawn as 
rectangles with one leg, and the combination of two legs into a 
line joining two tensors means contraction of the correspond- 
ing indices, (b) In the IMPS construction, the three-legged 
tensors are replaced by pieces of a cylinder, the vacuum state 
is represented by a half sphere, and the combination of two 
pieces corresponds to multiplication. The modes of the oper- 
ators Vs. {zi) correspond to vibrations of the cylinder pieces at 
different frequencies, and removing modes corresponds to re- 
ducing the number of frequencies at which the cylinder pieces 
can vibrate. 



mode has a special status and because we are particularly 
interested in the simplest situation (5 = and M = 1, 
which involves only one mode. 

Evaluating the expectation values, we find the explicit 
form 



iP{si,...,sn) 



l-5 + 55A\{zT-'A 

\n<m / 

E\ ^ L Zm 
SnSm > - 

^ q \ Zn 



, (5) 



where 6s is defined to be unity if X]^=i Sn = and 
zero otherwise. Note that the first factor in square 
brackets is the contribution from the zero mode, while 
the exponential is a product of M factors coming from 
modes 1 to M. In this and the next section, we choose 
Zn = exp(27rm/iV), which gives a translationally invari- 
ant state for 6 = 0. (Note that the definition in ([3]) does 
not impose \^p) to be translationally invariant in general.) 
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B. Entropy 

To investigate the properties of the state (O, we first 
compute the Renyi entropy S^^^ — — log[Tr(p|^)], where 
Pl is the reduced density operator obtained by tracing 
out spins L + 1 to N. A small calculation shows that 



,(2) 



E 

Sl,...,SiV s'-, 



E 



X V*(si,---,sl,Sl+i, 



>sjv)/V'*(si, 
>Sjv)/V'*(si, 



E 



Si, 



sn) 

-2 



(6) 



where the stars mean complex conjugation. The right 
hand side of this equation has the form of an expecta- 
tion value, and we compute it numerically by using the 
Metropolis Monte Carlo algorithm [isj . 

For an MPS, the entropy of a subchain of length L 
is limited by the logarithm of the virtual dimension, 
which is independent of L, and for critical systems the 
entropy grows approximately as ln(L) when the length 
of the subchain is short compared to the length of the 
complete chain [l^. The results for the state (O given 
in Fig. [5] shows that the entropy increases significantly, 
when the vertex operators are truncated to a finite num- 
ber of modes. This is a bit surprising because we go from 
infinitely many modes representing infinite dimensional 
Hilbert spaces to a finite number of modes representing 
infinite dimensional Hilbert spaces and shows that a sin- 
gle bosonic mode is, in fact, sufficient to obtain states, 
which contain a significant amount of entropy. We also 
observe that qualitatively similar results are obtained if 
the zero mode is or is not included. 

In addition to the pictorial representation of the tran- 
sition from M = 1 to M — >■ oo in Fig. ^a) , it is inter- 
esting to investigate the approximate scaling behavior of 
S)^' with L as a function of M. Inspired by the critical 
case M — > oo, for which the entropy is expected to fol- 



low the relation S 



(2) 



ln(7Vsin(7rL/iV)/7r)/4 -|- constant 



p^ , we plot the entropy as a function of sin(7rL/7V) in 
Fig. [3] For each value of M we fit a power law of the form 



s'f'' — a(sin{TTL/N))'^ -I- c to the data, where a, 6, and c 



are fitting parameters, and we plot the exponent 6 as a 
function of M. The fits show that the power law is a good 
approximation when AI is small provided the subchain 
does not consist of only a few spins. For intermediate 
M's there are some discrepancies, and for M = 50 and 
M — >■ oo, the exponent is close to zero. This reflects the 
transition to a logarithmic scaling for A/ — >■ cxd. Specif- 
ically, we find that s'l^ = 0.246 ln(sin(7rL/7V)) -I- 2.39 
provides a good fit to the data for M — oo, which is 
in accordance with the expression given above. Figure |3] 
provides a schematic summary of the scaling behavior of 
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FIG. 2. (Color online) Renyi entropy S)^ of a subchain of 
length L for A'^ = 4000, q = 0.15, various values of M, 
and (5 = 1 (a) and 5 = (b) obtained from Monte Carlo 
simulations. The faint dots on both sides of the curves for 
M = 1, 2, 5, 10, 20, 50 indicate the uncertainty in the numer- 
ical estimates. These curves are computed by averaging the 
right hand side r of ((6|) over the outcome of 50 or more Monte 
Carlo trajectories with different initial conditions, and the 

dots show — In |^mean(r) ± •\/var(r)/2j , where mean(r) and 

var(r) are the sample mean and the sample variance of r, 
respectively. We note that we use the same Monte Carlo tra- 
jectories to estimate Sl for all values of L. 



the entropy for MPS, the critical IMPS studied in p^ . 
and the IMPS considered in the present article. 

The dependence of the entropy on a is exemplified in 
Fig.[5]for M = 1 and 5 = 0, where the entropy is observed 
to increase with a, and a very similar behavior is found 
for M = 1 and 5—1. We can hence use a to fine-tune 
the entropy to a desired value. In Appendices |X] and |B1 
we show that the entropy for M = 1 and 5 = is given 
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FIG. 3. The same data as in Fig. [2ja), but plotted as a 
function of sin(7rL/A'^). The black solid lines are fits of the 
form 5^^' = a[sin(7rL/A'^)]'' + c to the numerical data (gray 
crosses), where a, b, and c are fitting parameters. The 
inset shows the exponent 6 (crosses) as a function of M 
{b = 0.0093 for M s- oo), and the dotted line is the fit 
b = 0.967 exp(-0.246M). 



f 2) 

FIG. 5. Renyi entropy Sj^ of a subchain of length L for 
— 4000, AI = 1,5 = 0, and various values of a as indicated. 
The gray crosses are obtained from Monte Carlo simulations, 
and the black solid lines show the values predicted by the 
approximate relation in Eq. ((T)). 

C. Correlation function 
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FIG. 4. (Color online) Schematic view of the scaling of the en- 
tropy of a subchain with the length of the subchain for various 
models. For MPS the entropy is bounded by the logarithm of 
the virtual dimension, which is independent of L. For the crit- 
ical case, the entropy scales with ln[A^ sin(7rI//A*')/7r], which 
is on the boundary between constant behavior and a power 
law, and for the IMPS considered in the present paper, there 
is a transition from the critical case to a power law scaling of 
the entropy with sin(7rL/A^) as M decreases towards one. 



approximately by the relation 



4^^ --In 



/ det(Mi + Ms) 
det(A/i 4- M2 -I- Ms) 



(7) 



in the limit of large N, where Mi, M2, and A/3 are 
the eight-by-eight matrices defined in Eqs. (jB3l) . (jB4l) . 
and (IBSp . To obtain this result, we use that z„ — 
exp{2Trin/N) has practically the same value for several 
nearby values of n when N is large and that the sum 
of Sn for such a group of n- values approximately follow 
a Gaussian probability distribution when spin configura- 
tions are chosen at random. The expression (jT]) is also 
plotted in the figure, and we observe an excellent agree- 
ment with the numerical results. 



Another important characteristic is the behavior of the 
correlation function between two operators acting on dif- 
ferent spins. Specifically, we shall consider 



(a,(0)a,(fc)) 



(8) 



where az{k) is the third of the Pauli operators acting on 
the kth spin, (Tz(0) = <Tz{N), and the expectation value is 
computed with respect to the atomic state ([5]) . Note that 
\ip{si, . . . , sn)\'^ is translationally invariant, and hence 
{(J z{n)a z{n + k)) = {az{0)az{k)), where n is an integer, 
and the argument of CTz is understood to be modulus N. 

Numerical results for the correlation function are given 
in Fig. |5] for various values of M , and we observe an 
excellent agreement with the approximate relation (for 
k N modulus N) obtained in App. |D] 



'5 + 2 Ei=i cos(2m7rfc/iV) 



{<Jz{0)<Jz{k)) = - 



N 



(9) 



which is expected to be valid for M ^ N. It is inter- 
esting to note that the correlation function for M = 1 
does not decay with the separation k for fixed N, but in- 
stead exhibits an oscillatory behavior with a wavelength, 
which is equal to the length of the chain. This is a very 
different behavior from the critical case A/ — > 00 and 
(5 = 1, in which the correlation function shows an anti- 
ferromagnetic ordering and a magnitude decreasing with 
k It is also in contrast to MPS, for which the corre- 
lation function decays exponentially with the separation 
[l| . We also note that the correlation functions for 5 = 1 
and 6 = only differ by the constant term —5/N, and 
both correlation functions do not depend on a. 
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FIG. 6. The correlation function (|8]) as a function of the sep- 
aration k between the two spins for A'^ = 500, a = 0.15, and 
S — 1 (a.) and S — (b). The gray crosses show results ob- 
tained from Monte Carlo simulations, and the solid black lines 
are the approximation in Eq. ((9)1. (Note that the point fc = 0, 
for which the correlation function is unity, is not included in 
the plots.) 
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FIG. 7. (Color online) To prepare the state ((5} experimen- 
tally, we initialize spins, which are here encoded in atoms, 
in the state (|-|-1) + | — 1))®^, and for each mode included in 
the vertex operators we let the spins interact one by one with 
a cavity field, which is subsequently subjected to a measure- 
ment. The inset shows the level structure of the atoms and 
the interaction with the cavity field and a classical field as 
described further in the text. 



imply 



N 



1^) = ( n 



(n)Va0Og'5o-j(n)VQTro ln(2„) 



n n 



\si,...,sn), (10) 



where 0^=i ^'-'^ some set of operators o„ is understood 
to mean 01O2 ■ ■ ■ on , and it is hence possible to prepare 
the desired state from a simple initial state if one can 
implement the M + S operators in ((T0|) . 

Let us first consider the operator coming from the mth 
mode. This operator is proportional to 



n 



exp 



,-1/2 



(11) 



III. EXPERIMENTAL IMPLEMENTATION 



Method 



The considered class of states is also interesting from 
the point of view that the structure of Eq. ^ and the fact 
that Vs„(2;„) only depends on the state of the nth spin 
provide a direct recipe for preparing the states experi- 
mentally in which an ancilla system sequentially interacts 
once with each of the spins. The form of the truncated 
vertex operators furthermore allows us to decompose the 
ancilla system into M -\- 5 independent bosonic modes 
that can be treated separately. Specifically, ([3]) and ^ 



which for z„ — exp(27rm/A^) is the expectation value of 
a product of N displacement operators, where the dis- 
placement az{n)^/am~^/'^ exp{2TTinm/N) accomplished 
by the nth displacement operator depends on the state of 
the nth spin. In other words, we can implement the con- 
tribution from the mth mode by initializing the bosonic 
mode in the vacuum state, engineering the Hamiltonian 



H = ihaz{n){gal^ - 9* am), 



(12) 



where g is a coupling strength, applying it for a time t 
such that gdt — -yam-^/^ exp(27rm7n/A^) for each of 
the spins, and finally projecting the bosonic mode on the 
vacuum state via a conditional measurement. 

The bosonic mode could be a single mode of a cavity, 
and the spins could be atoms traversing the cavity one 
by one as illustrated in Fig. [T] The movements of the 
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atoms in and out of the cavity could, for instance, be 
accomplished by the method demonstrated in [l^ , where 
laser fields are used to form a conveyor belt for the atoms. 
The two levels of the spins correspond to two ground 
state levels |±1) of the atoms, and we assume that the 
cavity field couples only the state | + to an excited state. 
By choosing the cavity field to be far detuned from the 
atomic transition and driving the same transition with 
a classical field at the same frequency, one achieves the 
Hamiltonian 

H = -^mf + ~gnal+rn*a„, + \~g\\la,r^)\ + l){+l\, 

,(13) 

where A is the detuning between the atomic transition 
and the frequency of the cavity mode, g is the coupling 
strength between the cavity mode and the atom, and il 
is the Rabi frequency for the driving of the atom with 
the classical field. The ac stark shift -fi,|f7p|+l)(+l|/A 
can be compensated by coupling | + 1) to another ex- 
cited state with another classical field, and the term 
— ?i|^paJjjam|-|-l)(-|-l|/A can be neglected if f2 is sufH- 
ciently large. The resulting interaction only displaces the 
amplitude of the cavity field if the atom is in the state 
I -1-1). One way to correct for this is to apply a displace- 
ment operation to the cavity field, which is independent 
of the state of the atom. This is, in fact, not necessary, 
however, because the error in the displacement occurring 
when the nth atom interacts with the cavity field is op- 
posite to the error in the displacement occurring when 
the (n-f iV/2)th atom interacts with the cavity field. Ex- 
cept for an overall phase factor, the only consequence of 
not applying the additional displacement operations is to 
introduce some single atom phase factors, but these can 
be canceled by applying the unitary operator 



N 



exp 2ia [<^z (n) + 1] Im 



^2ni{n-m)/N 
Km=l J 



(14) 



to the atoms after all the interactions with the cavity field 
have taken place. The final projection of the cavity field 
on the vacuum state may, for instance, be accomplished 
with a photo detector, or one can use atoms to probe the 
state of the cavity field (see [l3| for an implementation 
in the microwave regime). 

The contribution from the zero mode (if it is included) 
can be implemented along similar lines. The key point is 
to note that 



N 



lim ( T~r exp \/~^s. 



(a\i 



,a-nin/{l3N) 



floe 



-a-iTin/{liN) 



(15) 



where ao is a bosonic annihilation operator and the ex- 
pectation value is evaluated in the vacuum state of that 
operator. We can hence use the same procedure as be- 
fore if we take the limit where \^ g^ldt\ goes to infinity 



while the deviation of the phase of / gVLdt from — 7r/2 
goes to zero as 1/|/ g{ldt\^. The physical reason why a 
5s appears in this limit is that the norm of the sum of 
the displacements is either zero or infinite when the am- 
plitude of the classical field is infinite, and there is only a 
nonvanishing overlap with the vacuum state of the cavity 
mode after the interactions in the former case. 

We finally need a way to initialize the atoms in the 
state s kii ■ • ■ 5 sjv)- This can be done by pump- 

ing all the atoms to the state | — 1) and applying a 7r/2 
pulse between the levels | — 1) and |-|-1). 



B. Success probability 

The proposed generation scheme is only efficient if the 
probability Pm for successfully projecting all the cavity 
modes on the vacuum state after the interactions with 
the atoms is not too small. Pm is the square of the norm 
of the final state of the atoms relative to the square of 
the norm of the initial state of the atoms, and we find 



P, 



M 



where 



si,...,sjv \ g=l / 



1 ^ 



/N 



(16) 



(17) 



We evaluate the expression in ()16|) numerically by in- 
cluding a larger and larger number of randomly chosen 
configurations of the spins in the sum until the sum di- 
vided by the number of included configurations converges 
to a constant value, and the results are shown in Fig. [8l 
It is possible to obtain an analytical approximation to 
Pm as follows. In App. [X] we show that the values of 
= A^Re(si) and = A^Im(si) for randomly cho- 
sen spin configurations in the limit of large follow a 
Gaussian probability distribution 

P{xN,yN) = ^cxp[-{x% + y%)/N], (18) 

and the same distribution also applies to the real and 
imaginary parts of Nsq. Hence, if we start from the 
atomic state |si, . . . , sn) and only apply the op- 

erations required to take the qth mode in the truncated 
vertex operators into account, the projection on the vac- 
uum state for that mode proceeds successfully with prob- 
ability 



dx 



dVN^ exp 



xP{x'-^\y^^'^)=[l + 



(19) 
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FIG. 8. Probability of successfully projecting all the Al + 5 
cavity modes on the vacuum state after the interactions with 
the atoms as a function of the number of atoms for a = 0.15 
and 5 — 1 (a.) and 5 = (b). The gray crosses are numerical 
results, and the black solid lines represent the approximate 
relation in Eq. ((20)) . 



Since the inclusion of the zero mode projects the atomic 
state on the subspace of states with J2n=i = in 
addition to introducing some phase factors, the proba- 
bility of success when we only include the zero mode is 
2-^ X iV!/[(A^/2)!]2, which by use of Stirling's approxi- 
mation can be replaced by [2/ {nN)]^^"^ for large iV. 

From the arguments in App.[C]we expect the distribu- 
tions of Nsq for a moderate number of different values of 
q to be approximately independent in the limit of large 
N. This means that the success probability Pm factor- 
izes into a product of the individual contributions from 
the M -\- S modes, and we hence predict the approximate 
relation 



(20) 




9=1 



cellent agreement with the numerical computations even 
for only moderately high N. Equation (^0)) suggests that 
Pm scales approximately as iV^(^^+'^/^' , and for moder- 
ate M we hence conclude that the number of operations 
required to prepare the desired state grows only polyno- 
mially in A^. The proposed method can, however, not 
be used for generating the critical state investigated in 
[lo| since the success probability vanishes in the limit 
M -> oo. 



IV. MODIFICATIONS 

The experimental preparation scheme outlined above 
can be modified in various ways, which lead to new states 
that extend the family of states considered above, and we 
investigate the effects of two such modifications in the 
following. We focus on the case of a single bosonic mode, 
i.e., M = 1 and 6 — 0, throughout. 



Squeezing of the bosonic mode before and after 
the interactions 



One way to modify the preparation scheme is to modify 
the state of the cavity field before and after the interac- 
tions with the atoms. A nontrivial example is to apply 
a squeezing operation, i.e., the time evolution operator 
exp(r(e^*'^(a|)^ — e~^*'^ai)/2), where r and (f> are real pa- 
rameters, before the interactions and to apply the inverse 
operation after the interactions. Physically, this can be 
done by placing a nonlinear crystal inside the cavity and 
pumping it with a classical field (isj , and mathematically 
this amounts to replacing the vacuum expectation value 
in Eq. ([3]) by an expectation value with respect to the 
squeezed vacuum state 



|sq) - (l-4|Ani/V('^I)^|0), 
where A = e^^"^ tanh(r)/2, i.e.. 



(21) 



N 



(2-iTin/N) 



X e 



n=l 

— Sn^/aai cxp(— 27rm/A^) 



exp 



E 



a > 

n < rn 

21412 2aN^ 



sq) 

2-Ki(m-n) / N 



4a7V^ ,2 
1-4|A|2 l^^il 



l-A\A\- 



(22) 



This result is also plotted in Fig. |S1 and we see an ex- 



where si is defined in Eq. PT)) and we have used a theo- 
rem derived in [l^ to evaluate the expectation value. 

We can use the same approach as for A = Q to de- 
rive analytical approximations to the Renyi entropy (see 
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FIG. 9. Renyi entropy S^^^ as a function of the length of 
the subchain for A'^ = 4000, a — 0.15, and various values 
of squeezing of the cavity field. The gray crosses are results 
from Monte Carlo simulations, and the black solid lines are 
the analytical approximation derived in Add. [B1 



FIG. 10. The correlation function ((8]) as a function of the 
separation between the spins for = 500, a = 0.15, and 
various values of the squeezing of the cavity field. The gray 
crosses are results from Monte Carlo simulations, and the 
black solid lines are the analytical approximation in Eq. (|23|) . 



App.[B]), the correlation function (see App. [D|) 

{a^{0)a,{k)) = -cos(27rfc/A^) 

2a(l + 4|Ap)/(l - 4|Ap) + 2a^jV 
^ 1 + 2aA^(l + 4|A|2)/(1 - 4|A|2) + a^N^ 

(for k N modulus N), and the success probability 



(23) 



1 

2^ 



E 



exp 



aN 



1 ~4|A|2 



l + 4|Ap 
1-4|A|2 

Re {A*Sf) 



iSil^ 



1 + 2aN 



1 + 4|A|2 
1 -4|A|2 



-1/2 



(24) 



Numerical and analytical results for the entropy and the 
correlation function are provided in Figs. IHl and I10[ and 
it is apparent that \A\ needs to be close to 1/2 in order to 
observe significant deviations from the case A = 0. More 
precisely, Eqs. ^ and ^ show that (1 + 4|A|2)/(1 - 
4|A|2) should be on the order of aN or larger for the 
squeezing to have a significant effect. If \A\ is close to 
1/2, we observe an increase in the entropy compared to 
the case A — i). The correlation function keeps the same 
shape but the amplitude decreases from 2/N to \/N as 
A increases from to 1/2. Finally, we note that Pgq ~ 
^/7/JaN) when \A\^l/2-e and e < l/{aN). 



B. Projection of a part of the bosonic mode on 
vacuum between each interaction 



Another natural way to modify the preparation scheme 
is to allow a bit of the cavity field to leak out between each 



interaction with one atom and project the field that has 
escaped from the cavity on the vacuum state to remain 
within the set of pure states. The leakage corresponds to 
the application of a beam splitter operation, and hence 



N 



fcj,-gajb„ 



. n=l 



X e 



(25) 



where &„ are bosonic annihilation operators, i.e., 
\bn,h\j^ — Snm, and g is a parameter quantifying the 
amount of leakage (the relative loss in intensity is 
sin^d^D). Commuting operators in this expression, one 
can show that it is again possible to only displace the 
cavity field each time an atom is in the state |+1) pro- 
vided (i) one applies appropriate displacement operations 
to the field that leaks out of the cavity, (ii) one applies 
an appropriate displacement to the cavity field after com- 
pleting the interaction with all the atoms and before the 
projection on the vacuum state, and (iii) one applies ap- 
propriate single atom phase shifts to the atoms. We note 
that the displacements of the field amplitudes can be ac- 
complished by use of a strong classical field and a beam 
splitter as proposed in [2^. 

Evaluating the expectation value, we find that the 
wave function 



V'pr(si,---,SAr) 



n exp{-as„s™[cos(|ff|)]'"-"e2-("-")/^} (26) 

is, in fact, again given by Eq. ([SJ (with M = 1 and 
S = 0), but now Zn = cos(|(7|)"e2'^™/^. The numerical 
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FIG. 11. (Color online) Renyi entropy S^^' obtained from 
Monte Carlo simulations with A'^ = 4000, a = 0.15, and vari- 
ous amounts of leakage. 
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FIG. 12. (Color online) The correlation function {a-z{0)az{k)} 
obtained from Monte Carlo simulations with TV = 4000, a — 
0.15, and various amounts of leakage. 



results for the entropy and the correlation function ([8]) in 
Figs. fTD and [T2] show that even a small leakage has a sig- 
nificant impact on the properties of the produced states. 
The entropy is reduced significantly, and the correlation 
function now decays with separation. For g = 0.05 x tt/2, 
the decay is roughly exponential. This indicates that it 
is important to exclude losses from the setup to obtain 
the high entropies and long range correlations found in 
Sec. im The correlation function also shows that an- 
tialignment of nearby spins become more favorable when 
leakage is included. This can be understood from the fact 
that the field leaking out of the cavity is more likely to be 
projected on the vacuum state if the total displacement 
of the cavity field is close to zero at all times during the 
preparation. 



FIG. 13. (Color online) Success probability as a function of 
the number of spins for a = 0.15 and various values of g. The 
solid line almost following the points for g — 0.05 x 7r/2 is the 
fit Ppr = exp(-0.0704iV°-*227 _ o.6043). 

The probability of success 

1 r N N 

^P-- = 2^ j " " H 51 

Si,...,Sjsf \ 7l—lm — l 

X [cos(5)]l™-"l cos[27r(m - n)/iV]| (27) 

is plotted as a function of the number of spins in Fig. [131 
Ppr decreases significantly with increasing \g\, and for 
g = 0.05 X 7r/2, the decrease in Pp^ with N is almost 
exponential as detailed in the figure. 

V. CONCLUSION 

In conclusion, we have proposed and investigated the 
properties of a specific class of states of a one-dimensional 
spin chain, and we have suggested a (probabilistic) 
scheme for preparing the states experimentally, for which 
the required number of operations scales only polynomi- 
ally with the total number of spins N. The class of states 
is obtained by replacing the matrices in a matrix prod- 
uct state by chiral vertex operators of a bosonic confor- 
mal field theory and truncating these operators to involve 
only a finite number of bosonic modes. This construction 
is interesting, because it makes it apparent how one can 
generate the states by sequential interaction of the spins 
with an ancilla system, and because the resulting states 
turn out to have properties that are very different from 
MPS and also from the critical model obtained without 
truncating the vertex operators. 

Specifically, we have found that the Renyi entropy of 
a subchain of length L increases significantly, when the 
vertex operators are truncated, and the functional de- 
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pendence of the entropy on L turns out to approximately 
follow a power law in sin(7rL/7V) for small M as opposed 
to the logarithmic behavior observed in the critical case. 
We have thus demonstrated that replacing the finite di- 
mensional matrices in MPS with infinite dimensional op- 
erators acting on a single bosonic mode is sufficient to 
violate the area law and obtain states with high entropy. 
Considering the Uz-Uz correlation function for two spins 
separated by a distance fc, we have found undamped os- 
cillatory behavior as a function of k when keeping only a 
few modes, which is in contrast to the power law decay 
with k for critical systems and the exponential decay for 
MPS. These differences mean that the considered class 
of states provides a framework for describing a particu- 
lar set of problems, which can not be handled by MPS or 
critical models. Finally, we have demonstrated how the 
properties of the states can be altered by modifying the 
preparation scheme in various ways. 
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Appendix A: Distribution of X]^=i ^-nC' 



In this appendix, we determine the probability distri- 
bution in the limit of large N for obtaining given values 
of 



XL 



VL 



27Tin/N 



Im (^E^^e^-/^^ 



(Al) 
(A2) 



if we choose random configurations of the spins. For N 
large, we can to a good approximation replace the an- 
gle 4> = 2Trn/N by a continuous parameter and assume 
that several terms in each of the sums have a value of 4> 
within an interval of practically infinitesimal width d(j), 
i.e., dN = Nd(j)/ {2Tr) is large compared to one. Consid- 
ering one such group of dN terms with G (jjj + d(f>[, 
= (j — l)d0, the probability to obtain a given value 
of Wj = Srij , where rij labels the terms within the 
group, is 



dm 



2'"^ [{dN + Wj)/2]l[{dN ~ Wj)/2] 



(A3) 



which in the limit of large dN reduces to a Gaussian 
distribution 



with mean value and variance dN. Rewriting (IA1[) and 
(|A2p into XL = cos((/)j) and j/l = J2j 

we note that xl and j/l are linear combinations of ran- 
dom Gaussian variables, and xl and i/l are hence also 
random variables with a Gaussian probability distribu- 
tion P{xL,yL)- Since {xl) = '^j{'Wj)'^os{4>j) = and 
{ul) = siu((/)j) = 0, we conclude that 



P{xL,yL) 



27r^det(F) 



exp 



(A5) 

in the limit of large iV, where V is the covariance matrix 



V 



\xl) {xlVl) 
{xlVl) {vD 



(A6) 



L N 



cos\(t))d(t) = - + — sin Att— 



2 Stt 



N 



(A7) 



{yl) 



(xlVl) 



f^-^L/N L N f 

2^1 sin^(^W = 2-^sin(^4. 



N 
2^ 



27rL/N 



sm{(j)) cos((?!))c?0 



N 
8^ 



1 — COS 47r 



L 

N 



(A8) 



(A9) 



and we have used (wjWji) = for j ^ /. As a 
final remark, we note that the probability distribu- 
tion PiiL,PL) of S:l = Re (^^^^^^^ s^e^W^v) and 

PL = Im (^n=L+i SriG^'^™/^^ is a Gaussian with (xl) — 

(Pl) — and covariance matrix V with entries (xj^) = 
N/2-{xl), iPl) = N/2-{pI), and {xlPl) = -{xlPl), 
whereas (xlXl), (xlPl), (plxl), and {plPl) are aU zero. 



Appendix B: Analytical expression for the Renyi 
entropy for M — 1 and 5 = 

Combining ^ and © for M = 1 and i5 = 0, we obtain 



(xL+iyL~x'^-iv'L)(S:L-iyL-S:'^+iyi) 



X e 



E 



-a[(xL+iLf + (pL+PLf] 



, (Bl) 



1 



V2^KdN 



exp 



2dN 



(A4) where Xl, Pl, xl and pL are defined in App. [X\ and x^, 
p^, x'j^ and p'j^ are the same quantities evaluated for the 
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spin configuration s[, . . . , s^. For large N, the results in 
App. [Xl allow us to make the replacement 



where a. 



is the third of the Pauli matrices. 



1 

2N 



E 

Si,...,Sk 



dxj 



dxj 



/CO 
dyLP{xL,yL)P{iL,yL) (B2) 

and the same for the sum over Si, . . . , s'j^. The right hand 
side of (jBl[) then turns into a ratio between two Gaussian 
integrals that are easily evaluated to give Eq. ([7]) with 



Ml 



h h 

h h 

0/2/2 

0/2/2 



Appendix C: Approximate independence of 

E^^i s„e2-'"f/^ and E^Li Sne^"*"'/^ for p / g 

Assuming we can drag the Gaussian approximation in 
App. El to the limit, where each individual spin is re- 
placed by a Gaussian random variable Wn with mean 
and variance 1 , we can write the real and imaginary parts 
of J2n=i StiG^'^™''^^, where p is a positive 



and 



mteger, as 



(B3) 



Mo 



a 
'2 






7712 





-1712 


T 

m2 





T 

-1712 








-TO2 





7712 


T 

-1^2 





T 






Jp) 



N 



Wn cos(2TTnp/ N) , 



n=l 
N 



(B4) 



Pn' = ^WnSm{2TTnp/N), 



(CI) 
(C2) 



and 



and it follows that 



M3 = - 





























(B5) 



2^ 7o 



cos(p(/)) cos{q<j))d<j) — SpqN/2, (C3) 



2Tr 



and cTj, = 







— / sin(p0) sin(g0)d0 = (5pgAf/2, (C4) 

— / cos{p(f>) sm{q(j))d(j) = 0. (C5) 
27r Jo 

is the second of the Pauli matrices. por a Gaussian probability distribution, this is sufHcient 



where I2 is the two-by-two identity matrix, 7712 — 12+ fJ. 



Repeating the same arguments for the case of nonzero 
squeezing considered in Sec. IIV Al we again arrive at 
Eq. (O, but the matrices Mi and M2 are replaced by 



(p) 

N 



and p^'' are independent of 



Ml = a 



1 + 4|A|^ 
l-A\A\^ 



I2 h 

/2 /2 

0/2/2 

0/2/2 



4\A\ 



1-4|A|2 



fJ 
cr 








(B6) 



and 



M2 = 



"0 7712 - 


-7712 




Tnl' —771^ 







-7712 


7n,2 




—mj 7^1" 









" I2 





-/21 


4a|Ap 


h 


-h 





1 -4|A|2 


-h 





h 




. -h 


h 







" a. 







2a\A\ 


0, 


-<yz 





1-4|^|2 


-a. 









.-0, 








to conclude that x 

and pj^^ for p ^ q. We hence expect X]n=i Sne^'^"'^'/^ 
and X]n=i SnC^^™''/^ to be approximately independent 
for p ^ q. 

Another way to think of the independence of 
E^=i SnC^'^^f/^ and Y.n=i SnC^'^"''/^ is to regard the 
terms in the sums as points positioned at e'^^^^p/^ 
(^^2iTmq/N \^ in the complex plane and associated with ei- 
ther a minus or a plus depending on the value of s„. 
One can divide these points into groups with similar val- 
ues of 2'imp/N {2Tmq/N) modulus 27r, and the value 

of En=i Sne^'^^P/^ {T.n=i SnG^'^"'?/^) is determined by 
the number of points associated with plus and with minus 
within each group. By changing the phase of the points 
from 2-Knp/N to 2'imq/N , a reorganization of the division 
of the points into groups occurs, which depends on the 
precise value of 2'iTnp/N for the points associated with 
plus or minus belonging to the same group before the 
change. When there are many points within each group, 
it seems plausible that the values of Y^^^i Sne^'^™'^/^ at- 
tainable for a fixed value of '^2^=1 Sne^'^"'^/^ are almost 

except for 



independent of the value of X)^=i ■^"^ 



^2'Kinp / N 



very special choices of Y^^^i Sne^'^^"'^^^ , which justifies 
the assumption of independence. If, however, we consider 
several sums of the form X]^=i SnS^^^"''^^^ with different 
values of p, more constraints are present for the possible 



(B7) 
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values of one of them for all the others fixed. Hence, we 
expect the assumption of independence to break down 
if too many sums of the form J2n=i SnC^'^*"^/^ are in- 
volved. 



{az{0)crz{k)) correlation function. We first consider 
the case without squeezing. Since {o'z{0)o'z{k)) = 
{az{n)az{n + fc)), n ~ 1,2, ... ,N, it follows from ^ 
and ^ that 



Appendix D: Analytical expression for the 
((Tz (0)0-2 (fc)) correlation function 

In this appendix, we provide analytical arguments 
for the approximate expressions © and (1^51) for the 



{az{0)azik)) = 





-a J2^=i SnS„+m E^i ^ cos(2qnm/N) 


E.„.....„(l-'5 + <5<5.)exp 


E,^=l Em=l SnS„+m Egll ^ COs(2q7rTO/iV) 


N 



(Dl) 



This expression can be simplified by noting that 



N 



N 



q=l 



\sn\ 



N/2-1 

(-ir|Sjv/2|' + 2 ^ cos(27rm(z/A^)|S, 
9=1 



9I 



, (D2) 



N N ^ \ ^^1 

SnSn+m ^ - cos{2TTqm/N) = iV^ ^ -|s 



9I ' 



9=1 



9=1 



(D3) 



and (for k ^ N modulus N) 

N/2-1 

1 + 2 ^ cos{2T:kp/N) + (-1)*^ = 0, (D4) 
p=i 

where Sq is the discrete Fourier transform of s„ defined 
in Eq. i^Ji). For k ^ N modulus TV, ((DT|) reduces to 
(ctz(O)ctz(A^)) = 1, and for all other values of k, 

{aziO)azik)) = FMiN) - Fm{N/2) 
N/2-1 

-2 cosi2TTkp/N)[FM{N/2)^FM{p)], (D5) 
p=i 

where 

Eco„fk~Ppexp(-a7V2E:iiik^9p) 
Eco„fexp(-«iV2^,'lii|S,P 

is the expectation value of |spp and X^cont ^^'^ '^'^^'^ 
all allowed configurations of the spins. Note, in particu- 
lar, that the constraint X]^^=i ■^n = for 5 = 1 translates 



into sn = 0, whereas 

N 



n—1 9— 1 



1. 



(D7) 



If we assume Sp and Sq to be independent as discussed 
in App. [Cj we can write 

, V '^"-'^T^'^^.y for l<p<M 



for M <p< N/2 a.ndp = N 



The sum 



conf n—1 m—1 



1 27r'i(m-n)p/iV 



conf 



(D8) 

SnSm (D9) 



can be calculated exactly. For the case, where the zero 
mode is included, the total number of spin configurations 
is Nl/[{N/2)l]'^, the number of spin configurations with 
s„s„ = -1-1 is 2(iV-2)!/[(iV/2-2)!(Ar/2)!] for n ^ m, the 
number of spin configurations with s„Sm = —1 is 2 (TV — 
2)!/[(iV/2 - l)!(iV/2 - 1)!] for n ^ m, and consequently 



conf 



Iv:^ ((Af/2)!)^ forriT^m 



(DIO) 



For the case, where the zero mode is not included, the 
total number of spin configurations is 2^, the number of 
spin configurations with s„s,„ = -|-1 is 2^^^ for n 7^ to, 
the number of spin configurations with s„s„i = — 1 is 
2^~^ for n m, and 

Ef 2^ for n = m , . 

forn^m ■ ^^^^^ 

conf 

When we insert this into (|D9p . we obtain 



y fi 

Z— /conr 



l/{N-6) ioTp=l,2,...,N-l 
{1~S)/N ioTp = N 

(D12) 
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Note that this result is consistent with (ID7[) . 

Returning to the expression (|D8p for Fm{p) for 1 < 
we note that the exponential is small unless \sp\ 
is of order p/{aN'^) or smaller. The value of Fm{p) for 
1 < p < M is hence at least about a factor of N smaller 
than the value of Fm{p) for M < p < N/2. In conclusion, 
we thus have 

( for 1 < p < M 

Fm{p) ~ I for A/ + 1 < p < N/2 , (D13) 

[ {l-6)/N ioi-p = N 

which, when combined with (jD5[) . leads to (0). We note 
that (jD13|) is, in fact, inconsistent with (ID7p . and the 
inconsistency grows with M . This indicates that the ap- 
proximation is less accurate for higher M as expected. 

For the case of nonzero squeezing investigated in 
Sec. IIV Al the derivation follows similar steps. In par- 



ticular, Fi{p) is still given by (ID13I) for p 1. For p = 1, 
however, we can no longer assume that i^i(l), which in 
this case is given by 

^i(l) = 

E.i, exp (-ai^a;2^ - a^j/^y%^ 

(D14) 

is approximately zero {xn and are defined in (|Aip 
and (jA2p . respectively). In fact, for \A\ (l/2)~, only 
configurations with yN = contribute, and for these con- 
figurations the exponential factor is unity, while xn may 
take different values. Instead, we evaluate the expecta- 
tion value by utilizing the probability distribution for x m 
and pn derived in App. \^ and this leads to 
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